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Abstract. We construct an explicit model where it can be established if a two mode pure Gaussian system 
is entangled or not by acting only on one of the parts that constitute the system. Measuring the dispersion 
in momentum and the time evolution of the dispersion in position of one particle we can tell if entanglement 
is present as well as the degree of entanglement of the system. 

PACS. 03.67.-a Quantum information - 03.65Ud Entanglement and quantum nonlocality 



] 1 Introduction 

, One of the most intriguing features of Quantum Mechan- 

■ ics (QM) is entanglement and in the early times of QM 
' it was recognized by Erwin Schrodinger and by Ein- 
. stein, Podolsky and Rosen [21. Later, John Bell |3] showed 

■ that the non-local aspect of entanglement is experimen- 
I tally testable through his famous inequality. 

In recent years the interest in entanglement has in- 
' creased considerably. First because it is a fundamental tool 
I in Quantum Information Theory and a consistent char- 
' acterization of its theoretical properties is needed. Sec- 
' ond because the present stage of technology permits us to 
1 perform some experimental manipulations with it such as 
Quantum Teleportation 0| and Quantum Cryptography 

m- 

In the study of the properties of entanglement Peres |Zj 
and the Horodecki family |H] have derived a necessary and 
sufficient condition for the separability of 2 x 2 and 2x3 
systems. Some years later Simon ini and Duan et al [El 
have obtained a necessary and sufficient condition for the 
separability of two-party Gaussian states. Given that the 
state is non-separable we should have a measure of the de- 
gree of this inseparability. There are at least three distinct 
measures of entanglement: the entanglement of formation 
|12| . the distillable entanglement ^S] and the relative en- 
tropy of entanglement |14I15| . In any calculations done 
with these three measures of entanglement and the two 
criteria for separability we must use the total density ma- 
trix of the bipartite system. That is, given the density 
matrix that describes the whole bipartite system we can 
determine if the system is separable or not and its degree 
of entanglement. 

For pure states it is well known that the knowledge 
of the whole reduced density matrix allows us to decide 



whether or not the system is entangled. If the reduced 
density matrix is pure (Ti(pi) — 1) the system is separable 
and entangled if it is not pure (Tr(pf) < 1). This article 
aims to show that we need not know the whole reduced 
density matrix pi (or equivalently P2) of a bipartite pure 
system to deduce if it is entangled or not. Here we show 
that only the diagonal elements of the reduced density 
matrix, even if it is written in a representation where it is 
not diagonal, are sufficient to detect entanglement. 



To explicitly demonstrate this we construct two paradig- 
matic cases: a non-entangled two particle Gaussian wave 
function in configuration space and an entangled two par- 
ticle Gaussian wave function. We then let both systems 
freely evolve in time. We show that when studying an in- 
dividual particle of each case we get different results for 
the time evolution of the dispersion in position. This fact 
allows us to tell if we are working with a non-entangled 
or an entangled bipartite Gaussian wave function, and in 
the case of an entangled system we can also extract from 
this evolution the degree of entanglement. 



Some aspects of this approach are similar to the to- 
tal wave function reconstruction shown in Refs. |l(ill7ll8] . 
Here, however, we do not need to reconstruct the whole 
wave function describing the two particles. We only need 
two elements of the reduced wave function describing a 
single particle, i. e., its dispersion in position and in mo- 
mentum. We should also mention that by employing the 
powerful techniques given in Refs. |l(ill7| it might be pos- 
sible to generalize the following approach to the case of 
mixed Gaussian states. 
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2 The non-entangled bipartite system 

Consider a normalized one-dimensional separable two par- 
ticle Gaussian wave function where we assume, with no 
loss of generality, that the two particles have the same 
mass m but can in principle be distinguished from each 
other: 



3 The entangled bipartite system 

Let us now construct a normalized one-dimensional en- 
tangled two particle Gaussian wave function where we 
assume again, with no loss of generality, that the two 
particles have the same mass m but can in principle be 
distinguished from each other. 



where 



(1) 



(2) 



(3) 



Here uj{k) = is the dispersion relation for a free par- 
ticle and f{ki), /(— ^2) represents the fact that we have 
Gaussian particles moving in opposite directions |19| : 



fik) 



(27r)3/4 



■-(k~k,Y 



(4) 



In Eq. a represents the dispersion of the Gaussian 
wave packet centered in and the factor that multiplies 
the exponential is the normalization constant. 

Integrating in ki and k2^ then multiplying ^{xi,X2,t) 
by its complex conjugate and finally integrating in X2 we 
get the probability density of particle 1 at time t [T^: 



exp 



2 (xi - v^tf 



a2 l + F{t) 



where F{t) and Vc are given by: 
Fit) ^ Fit, a) ^^-.^ 



(5) 



Vc = . (6) 

m 



With Eq. ||SJ) we can calculate the dispersion Axi = 



(a;i) of the position of particle 1: 



Ax^it)^-,/lTFit)- 



(7) 



We can also obtain the dispersion of the momentum of 
particle 1 if we take the Fourier transform of Eq. . Then 
multiplying the result by its complex conjugate and inte- 
grating in k2 we obtain: 



■ exp 



-(fci 



(8) 



Using Eq. (jS)) and the fact that pi — Hki we easily get: 



Apiit) = -. 



(9) 



As expected for a free particle the dispersion in momen- 
tum is constant in time. 



>I'ixi,X2,t ^ 0) = J dkidk2f{ki,k2) 

x-0i(a;i,O) «) V2(a;2,0). 
Here ipiixi, 0) and "02(2^2, 0) are given by: 



(10) 
(11) 

V'2(a;2,0) -e^'^'^^^e^. (12) 

Eq. Hl()|l is a superposition of bipartite Gaussian wave 
packets centered in fci and k2 where /(fci, k2) — giki, /C2) 
5(fci -|- ^2) are the expansion coefhcients and 5(fci + ^2) is 
a restriction which entangles the system. This delta func- 
tion can be viewed as the requirement for the conservation 
of momentum in the center of mass frame, that is, we su- 
perpose bipartite Gaussian wave packets where each party 
moves in opposite directions centered at the same momen- 
tum. Eqs. Hll|l and (|12(l are proportional to Eqs. ^ and 
(j3Jl where we integrate for t = and substitute kc by fci 
and k2 respectively. By using the delta function Eq. ((Tn|l 
can be rewritten as: 



tZ/(a;i,X2,0) 



dkigiki 



(13) 



Eq. H13|l clearly shows that (5(fci -I- ^2) entangles our sys- 
tem. Because ]Pixi, X2,0) cannot be written as a simple 
tensor product of a wave function belonging to particle 
1 and another belonging to particle 2 we now deal with 
a non-separable wave function. Only if giki) is another 
delta function we can disentangle the system and recover 
Eq. J^l- In Eq. ((T^ g{ki) is chosen to be a Gaussian dis- 
tribution centered in kr'. 



-(6/2)2(fci - kef 



(14) 



where 6 is a new parameter that measures the degree 

2 

of entanglement as explained below and /„ = 1 -I- rip-, 
n — 1,2. We can see that when b ^ oo the function 
^exp [-(6/2)'(fci - kef] ^ 5(fci-fc,) ism and /2 -> 1, 
showing that entanglement has disappeared. This can be 
seen doing a straightforward calculation using Eqs. ((TUl 
and ifT^ : 



lim 'Fixi, a;2, 0) 



1/4 



1/4 



(15) 
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Eq. H15|l is identical to Eq. if we calculate the integrals 
in Eqs. (0 and JSJ- Furthermore, it can be shown that if 
6 — > and a — *■ oo Eq. l(T^ is the EPR state with xq — 
pi] . As stated in ref. Eq. l(T^ can be viewed as a 
generalized version of the EPR wave function. These two 
facts suggest that b should be considered as a measure of 
the degree of entanglement, where 6 ^ oo represents no 
entanglement and b represents the maximally entan- 
gled state. 

Doing the integral in Eq. H13() we get the normalized 
bipartite Gaussian wave function at t = 0: 



f 4, ika{xi-X2) 

2 •'2 6 



X exp 



■^(Xi +X2) + -^XiX2 



(16) 



It is interesting to note that Eq. H16|l represents a non- 
separable (entangled) state due to the term exp [^a;ia;2] . 
If — > 00 this term tends to 1 and we obtain Eq. ^ as 
a limiting case of Eq. IjlGI) . In other words, when 6 — )■ 00 
we have Eq. a separable, non-entangled state, and for 
any other value of b we have Eq. (|16|l . a non-separable, 
entangled state. 

In order to make rigorous that b furnishes the degree 
of entanglement of the state given by Eq. (|16|l and that 
only when & — > 00 we have a disentangled system we first 
calculate its correlation matrix (CM) and apply the Simon 
separability criterion j^, which shows that the bipartite 
Gaussian system is separable iff 6 ^ 00. After applying 
the Simon criterion, we make a local symplectic transfor- 
mation in the CM to put it in its standard form '[)'10' 
and then calculate its entanglement of formation (EoF) 
\22\ . which is a monotonically decreasing function of the 
parameter b, proving that the higher b the less entangled 
is the state. As we deal with a pure state, we note that 
we can calculate the von Neumann entropy of the reduced 
density matrix to obtain the entanglement of this system. 
However, we prefer using the EoF as given in Ref. [22] 
since, in the particular case of symmetric bipartite Gaus- 
sian states, it is more straightforward than the usual pro- 
cedure for pure states. 

The CM completely specify a two mode Gaussian state 
and it is a 4 X 4 matrix, which has the following elements 

mni: 

7y = Tr [{R,Rj + RjR,)p] - 2Tr[i?,p]Tr[i?jp], (17) 

where R — (Xi, Pi, ^2, ^2)"^ and Rj arc the position and 
momentum operators of the two particles. Doing the cal- 
culations we get the following CM: 



A C 
A 



(18) 



The Simon separability criterion says that the above CM 
represents a non-entangled system iff |ll) : 



/ = det^detS + {h^ ~ |detC|)^ 
-h^{detA^detB) > 0, 



Ti{AJCJBJC^J} 



(20) 



where J 



^Yoj- But a simple calculation shows that 
the rhs of Eq. ^ is: 



/ = -4.<l. 



(21) 



Hence, / < except when 6 — > cxd, proving that for any 
other value of b we have an entangled state. 

We now make the following local symplectic transfor- 
mation S = diag (s, s~^, s, s~^), where s = (4/1^/2/0'')^^^. 
This brings 7 to its standard form 70 = S'^S'^ |9I10| : 



70 




(22) 



where n = S/i/v^ and = kp = ?ia^/(&^v^)- This is 
a symmetric Gaussian system and Giedke et al. ,22, have 
shown that the EoF for this state is: 



EoF{9) = f 



\J{n- kx){n - kp) 



where, 

f{5) = c+(5) \og^[c+{5)] - c^{5) log2[c_(<5)]. 



(23) 



(24) 



Here c±{5) = ((5~^/^ ± (5^/^)^/4. Analyzing the behavior 
of the EoF given by Eq. H23|) we clearly see that it is a 
decreasing function of the parameter b (Figs. ^ and |2Jl. 




where 



Fig. 1. Entanglement of formation, Eq. 123^ . for the sym- 
metric Gaussian state given by Eq. 1161 . as function of the 
parameters a and b, where we have put ft = 1. We clearly see 
the EoF increasing as 6 ^ and decreasing as a ^ 0. 



Working in the Heisenberg picture we easily obtain for 
a free evolution. 
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Fig. 2. Entanglement of formation, Eq. 12311 . as a function of 
1/6 for ten values of the parameter a. From bottom to top the 
parameter a varies from 1 to 10 in increments of one unit. We 
have set h = 1. We clearly see that the EoF increases as b 
decreases and that for a given b, the higher a, the greater the 
EoF. 



^Piit) ^-VTi- (26) 
a 

Again, due to the free evolution of particle 1 the dispersion 
in momentum does not vary in time. 

We should mention that the formal solution of the 
Heisenberg equations of motion for the observables xi{t), 
xl{t), pi{t), and pl{t) are identical for the entangled and 
non-entangled case. Only when we take the mean values 
{xi{t)), {xl{t)), {pi{t)), and {Pi{t)) we obtain different 
quantities. This is due to the fact that we have differ- 
ent initial wave functions. In other words, entanglement 
manifests itself furnishing different initial conditions for 
the Heisenberg equations of motion, which imply different 
evolutions for the dispersions. 



4 The measuring protocol 

As we have all tools now, that is, all the dispersions in 
position and in momentum for the entangled and non- 
entangled case, we develop a measurement procedure to 
be used in an ensemble of two particle Gaussian systems 
which allows us to locally decide whether or not the par- 
ticles are entangled. From now on we assume H ^ m = I 
for simplicity. 

Let Bob be our physicist who receives one of the par- 
ticles of the bipartite Gaussian system produced by Al- 
ice. Bob knows, because Alice has told him, that all the 
particles he receives are either entangled or non-entangled 
Gaussian wave packets, according to the two constructions 
explained above. There are no other possibilities. Alice 
produces many pairs at once. And continues to produce 
many pairs at once for different times. Of course Bob does 
not know the values of the parameters a and b used by Al- 
ice. But Bob is curious enough and wants to know whether 
his particles are entangled or not. Bob cannot use any fur- 
ther classical communication, he can act only locally on 
his particles and he is able only to measure the dispersions 
in position and in momentum of his particles, that is, the 
diagonal elements of the system reduced density matrix. 
He proceeds as follows: 



First he measures, using a sub-ensemble, the dispersion 
in momentum of his wave packets. He obtains Api — u. He 
does not know yet whether Eq. lO or Eq. (|26|l represents 
what he measures. However he knows that it must be one 
of these two possibilities, which imply only two possible 
time evolution for the dispersion in the position of his 
particles. 

If his particles are not entangled and Bob uses in Eq. 
(UI) the fact that Api{t) = u = ^ he gets: 

Axiit) = — ^/l + AuH^. (27) 
2u 

But if Bob's particles are entangled and now he uses the 
fact that Api{t) = u = Eq. (j^SJl becomes: 

= (28) 

Looking at Eqs. H27I) and H28|) we see that if Bob knows 
at what time Alice has begun to produce the pairs he 
is able to discover, with only one measurement of Axi^ 
whether his particles are entangled with Alice's or not. 
The reason for this is simple: Let us suppose, with no loss 
of generality, that Alice begins to produce the pairs of 
particles at t — Q. Measuring the dispersion in position 
for a given time t Bob obtains Axi (t) . Remembering that 
Bob also knows the value of u, he can calculate, using 
Eq. I|27() . the value of Axi{t). If this calculated value of 
the dispersion agrees with the measured one, Bob has the 
non-entangled case. If these values of Axi are different, 
Bob has entanglement. In this last case, using Eq. (|28|l 
Bob can obtain the parameter b. For any t Bob can use 
this procedure. Bob sees two distinct curves for the time 
evolution of Axi{t), whether his particles are entangled 
or not. See Fig. 1 below: 



ax 




1 2 3 4 ^ ' 

Fig. 3. The dashed curve is the time evolution of the dispersion 
in position for an entangled Gaussian wave packet while the 
solid curve represents the non-entangled case. We have chosen 
6 = 1 and u = 1.01. 

Analyzing Eq. (|28|l we see that for it to be valid for all 
i > we must have for the entangled case: 

ub > 1. (29) 

It is worth noting that asymptotically Eq. (|27|l and (|28|l 
are the same. Therefore, in order for Bob to correctly dis- 
tinguish between the two cases he should make his mea- 
surements for times smaller than a critical time tc, which 
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is defined to be the time where the time independent term 
inside the square root of Eq. H28() is of the order of the 
term: 

~ 2Vu%^ - 1 ■ ^^^^ 
We can increase tf. making ub 1. This might seem as a 
limitation of our procedure but as Alice sends a classical 
message to Bob defining the origin of time, Bob can start 
making measurements as early as possible. 

Now let us make things harder to Bob. We assume from 
now on that Bob does not know when and where Alice 
has begun to produce the pairs. This fact means that Bob 
cannot use the previous procedure to answer whether or 
not his particles are entangled with Alice's. The previous 
protocol fails because Bob does not know what time t he 
should use to calculate Axi{t), which would have allowed 
him to compare this calculated value with the measured 
Axi{t). 

We first prove why a single measurement at time t is 
not enough for Bob to tell whether his particles are en- 
tangled or not. We are now assuming that he does not 
know when Alice has begun to produce the particles. The 
proof is achieved showing that the diagonal elements of 
the reduced density matrix (in position and in momentum 
representation) of the non-entangled system can be made 
identical to the diagonal elements of the reduced density 
matrix of the entangled system for t = 0. (The same rea- 
soning applies to any t, but for t — the calculations are 
much simpler and we do not lose in generality). 

For t = the diagonal elements of the reduced density 
matrix of the entangled system in momentum representa- 
tion is 

ei(fcl) - J {h,k2\^){1'\ki,k2)dk2 



exp 



(31) 



For any t, the diagonal elements of the reduced density 
matrix of the non-entangled system in the momentum rep- 
resentation, according to Eq. (|SJ, reads: 



pi(fci,i) 



{ki,k2\il^){ip\ki,k2) dk2 



-^ih-k^f 



(32) 



If we want identical diagonal elements of the reduced den- 
sity matrices we must impose that: 

, a 1 



(33) 



The diagonal elements of the reduced density matrix 
for t = of the entangled system written in the position 
representation is: 



{xi,X2\'l'){'l'\xi,X2) dx2 



2/2 

7ra2/] 



exp 



^x' 



(34) 



As stated in Eq. Q , the diagonal elements of the reduced 
density matrix for any t of the non-entangled system in 
position representation is: 

Pi{xi) = / {xi,X2\tp){i'\xi,X2) dX2 



7ra'2 ^1 + F(i,a') 

2 {xi - v^tf 



X exp 



a'2 1 + F(t, a') 



(35) 



If we want Eqs. (|34|l and H35|l giving the same statistical 
predictions we must have: 



2/2 ^ 2 1 
a2/i a'^ 1 + F{t, a') ' 



(36) 



Eq. 1)36(1 is a restriction which forces the two density ma- 
trix to give the same dispersion in position. (We do not 
need to bother with the first order moment of these Gaus- 
sian functions because a translation of the xi-axis sets it 
to zero.) If we use Eqs. (|33f) and H36|) and the fact that 
/„ = 1 + np- we arrive at the following condition: 



t 



1 



1 



2u2 Vu464 - 1 



(37) 



Eq. ((S7jl says that for only, and only one time t, the diago- 
nal elements of the reduced density matrices, one obtained 
from the entangled system and the other one obtained 
from the non-entangled system, furnish the same statis- 
tical predictions. This implies that single measurements 
of the dispersion in momentum and in position of parti- 
cle 1 do not tell us unequivocally whether we are dealing 
with a non-entangled or entangled Gaussian bipartite sys- 
tem. (Unless, of course, we know when Alice has begun to 
produce the pairs.) See Fig.^ 

To circumvent the limitation of the previous protocol 
Bob may apply the following one, which explicitly uses 
the difference in time evolution of the two systems: 

Bob again initially measures the dispersion in momen- 
tum of his particles {Api = u). As he does not know when 
and where Alice begins to produce the pairs of Gaussian 
particles, the time evolution of the dispersions in position 
for the non-entangled and entangled systems are: 



Axi{t) 



Ax^it)^- 



1 

2u 



1%^ 



4u4(t + to)2. 



(38) 



(39) 



Here tg is the time elapsed from the production of the pair 
by Alice until Bob makes his first set of measurements. 
Bob now makes several measurements of the dispersion in 
position for different times t. With these measurements he 
obtains the following set of points: 

{(Z\xi(0),0), {Axi{h),h), {Axi{tn),tn)}. 
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Fig. 4. The dashed curve is the time evolution of the dispersion 
in position for an entangled Gaussian wave packet produced 1 
unity of time after the production of the non-entangled case, 
which is represented by the solid curve. The curves intercept 
each other for t ~ 2.46. If Bob measures Axi for this time, 
he cannot distinguish between the two ways in which Alice 
can produce the pairs of particles. For any other point of the 
dashed curve we can find a solid one that crosses it. Therefore, 
Bob cannot distinguish how his particles were produced if he 
measures Axi only once. Here we have chosen 6=1 and u = 
1.01. 

He makes as many measurements as possible. With the n 
pairs of points above he fits the following curve, where a 
and (3 are the free parameters and u is already known: 



^2^1 W ^ — Va + 4w''(t + /3)2. (40) 
2u 

Looking at Eqs. gOJ), and we see that if the pa- 
rameter a — 1 Bob is deahng with non-entangled Gaussian 
functions, but if a 7^ 1 Bob deals with entangled particles. 
And using a Bob can calculate the value of the degree of 
entanglement b. Just for completeness we mention that 
P furnishes the time tg. For this protocol to be optimal, 
Bob should begin his measurements as soon as possible 
since, asymptotically in time, Eqs. H38|) and H39|l are seen 
to become identical. 



5 Conclusion 

We have shown an explicit model using two particle Gaus- 
sian systems where we can decide if we are dealing with 
non-entangled or entangled pairs acting only on one of the 
particles and measuring only the diagonal elements of its 
reduced density matrix. Measuring the dispersion in mo- 
mentum and then the time evolution of the dispersion in 
position of one member of the pair it is possible to dis- 
cern between the entangled and non-entangled cases. It is 
also possible with this procedure to determine the degree 
of the entanglement of the system. The above model sug- 
gests that just one part of the whole system can furnish 
more information about the degree of the entanglement of 
the system than we had imagined. 

Finally it is important to note that the presented mea- 
surement protocol uses the time evolution of the diago- 
nal elements of the reduced density matrix to determine 
whether or not we have entanglement. This fact shows that 
we may have a new tool to analyze the properties of entan- 
gled systems, i. e., the dynamical evolution of entangled 



states. So far all the methods used to study the properties 
of entangled systems have not employed the dynamics of 
the system. We are hopeful that studying the dynamics 
of entangled systems will help us to deepen our under- 
standing of entanglement and possibly it will unravel new 
features of entanglement not yet explored. 

Acknowledgments 

This work was supported by Fundagao de Amparo a Pes- 
quisa do Estado de Sao Paulo (FAPESP) and Conselho 
Nacional de Desenvolvimento Cientffico e Tecnologico (C- 
NPq). We thank Dr. Lea F. dos Santos for her careful 
reading of the manuscript and useful discussions. 

References 

1. E. Schrodinger, Proc. Camb. Phil. Soc. 31, 555 (1935). 

2. A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev. 47, 
777 (1935). 

3. J. S. Bell, Physics 1, 195 (1964). 

4. C. H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, A. 
Peres, and W. K. Wootters, Phys. Rev. Lett. 70, 1895 
(1993). 

5. C. H. Bennett, G. Brassard, and N. D. Mermin, Phys. Rev. 
Lett. 68, 557 (1992). 

6. A. K. Ekert, Phys. Rev. Lett. 67, 661 (1991). 

7. A. Peres, Phys. Rev. Lett. 77, 1413 (1996). 

8. M. Horodecki, P. Horodecki, and R. Horodecki, Phys. Lett. 
A 223, 1 (1996). 

9. R. Simon, Phys. Rev. Lett. 84, 2726 (2000). 

10. L-M. Duan, G. Giedke, J. I. Cirac, and P. ZoUer, Phys. 
Rev. Lett. 84, 2722 (2000). 

11. To recover the Simon separability criterion as it was pre- 
sented in Simon's paper [21 we should substitute in 
Eq. 12011 by ft^/4. This difference is due to the fact that 
Simon's CM is defined as the half of our CM 7. We use 
the same notation as |10| . 

12. C. H. Bennett, D. P. DiVincenzo, J. A. Smolin, and W. K. 
Wootters, Phys. Rev. A 54, 3824 (1996). 

13. C. H. Bennett, G. Brassard, S. Popescu, B. Schumacher, 
J. A. Smohn, and W. K. Wootters, Phys. Rev. Lett. 76, 
722 (1996). 

14. V. Vedral, M. B. Plenio, M. A. Rippin, and P. L. Knight, 
Phys. Rev. Lett. 78, 2275 (1997). 

15. V. Vedral, M. B. Plenio, K. Jacobs, and P. L. Knight, Phys. 
Rev. A 56, 4452 (1997). 

16. U. Leonhardt and M. G. Raymer, Phys. Rev. Lett. 76, 
1985 (1996). 

17. T. Richter and A. Wiinsche, Phys. Rev. A 53, R1974 

(1996) . 

18. Ch. Kurtsiefer, T. Pfau, and J. Mlynek, Nature 386, 150 

(1997) . 

19. C. Cohen-Tannoudji, B. Din, and F. Laloe, Quantum Me- 
chanics, Vol. 2 (Hermann and John Wiley & Sons, Paris, 
1977), Chap. 1, Complement Gi. 

20. G. Arfken, Mathematical Methods for Physicists, (Aca- 
demic Press, Orlando, 1985), Chap 8. 

21. B-G. Englert and K. Wodkiewicz, P hys. Rev. A 65, 
054303 (2002). 

22. G. Giedke, M. M. Wolf, O. Kriiger, R. F. Werner, and J. 
I. Cirac, Phys. Rev. Lett. 91, 107901 (2003). 



